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Abstract: This paper proposes a hybrid soft tissue model that consists of a multilayer structure and many spheres for sur- 
gical simulation system based on meshless. To improve accuracy of the model, tension is added to the three-parameter 
viscoelastic structure that connects the two spheres. By using haptic device, the three-parameter viscoelastic model (TPM) 
produces accurate deformationand also has better stress-strain, stress relaxation and creep properties. Stress relaxation 
and creep formulas have been obtained by mathematical formula derivation. Comparing with the experimental results of 
the real pig liver which were reported by Evren et al. and Amy et al, the curve lines of stress-strain, stress relaxation and 
creep of TPM are close to the experimental data of the real liver. Simulated results show that TPM has better real-time, 
stability and accuracy. 
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1. INTRODUCTION 

Virtual surgery simulation enables physicians to interact 
with flexibly customized simulate environments based on 
visual, auditory and haptic feedback without potentially 
harmful contact with real patients. With the use of virtual 
surgery, surgeons have the opportunity to surgically train in 
a low-cost and ethically sound environment [1]. On the other 
hand, real-time surgical simulation and medical training 
must employ the physical modeling and behavioral modeling 
of soft tissue so as to achieve realistic visual and tactile ef- 
fects [2]. 

Deformable models can be classified into two categories: 
physics-based and geometry-based. Physics-based methods 
are based on continuum mechanics, and can get accurate 
simulation results by directly solving the partial differential 
equations (PDEs) using numerical methods. Some of the 
prevailing methods include the finite element method 
(FEM), mass-spring model (MSM), boundary element 
method (BEM), point-based method. Geometry-based mod- 
els use intuitive methods instead of solving PDE. For exam- 
ple, vertex based model and spline based model. The advan- 
tages of geometry-based approaches are: fast and smooth, 
but the major drawback is in their accuracy. In contrast, 
physics-model is more accurate than geometry based ap- 
proaches, but it has the disadvantages of being slow and 
rough. The key point to build a physical model of soft tissue 
is to reduce computational complexity and simultaneously 
keep the accuracy of deformation effects [2]. 
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A material is called viscoelasticity if the stress-strain 
relationship of the material depends on time under applied 
forces. Ratchada [3] applied the method of tissue creep and 
stress relaxation to parallel-fibered collagenous tissues. Co- 
tin [4] established a hybrid elastic model for surgery training 
and simulation. Suzuki [5] created a deformable organ model 
by using sphere-filled method instead of the finite element 
method. But the deformation of the model is still based on 
MSM. Recently, meshless methods have been introduced 
in the computer graphics field for a real-time simulation of 
deformable objects. The meshless method, in contrast to the 
mesh-based method, discretized the object by merely using a 
set of nodes without a mesh structure. This method thus 
avoided the problems associated with remeshing and handles 
the topology change of the object in a relatively simple man- 
ner. The topological relation among meshless simulation 
nodes is defined by support of the node. Adjacent nodes, 
which lie within the support, acted like the nodes which are 
connected by the edges of the mesh in the mesh-based 
method [6]. Then researchers start to set up the meshless 
physical model such as potential energy method [7], discrete 
mechanics [8], which has become a popular research in soft 
tissue modeling and cutting. 

Haptic feedback is especially important in delicate surgi- 
cal pressures requiring fine sensations, especially when 
slightly excessive pressure can injure a patient. The ability to 
touch the virtual objects in the environments is one factor 
that provides a more realistic experience and is particularly 
suited to virtual training applications [9]. So, the create proc- 
ess is related to the constraint force needed to keep a valid 
model which is close to the real tissue from the relaxation 
and creep and other properties. 
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2. MODEL STRUCTURE 

2.1. Layered Sphere-Filled Structure 

A multilayer structure has been adopted to arrange the 
TPM, as shown in Fig. (1, a). The TPM has 424 spheres and 
the radius of the sphere is 10mm. These spheres are con- 
nected by a three-parameter structure which is composed of a 
Kelvin model and a linear spring. 

Under external force, the deformation of the outer-layer 
sphere is obvious, and the deformation of the inner-layer 
sphere will weaken orderly. The deformation effects fit the 
deformation of the real soft tissue, which is shown in Fig. (1, 
b). Fig. (2, a) presents the structure of layered sphere-filled 
based on meshless method in TPM. This model has been 
established by using VC++ and the OpenGL. 

To clarify this point further, the liver tissue model has 
been established by using MSM structure in this paper. 
Loading 2 newtons (N) in the Z-axis positive direction which 
act on the MSM, the deformation effects of the liver model 
are shown in Fig. (3, a). Similarly, loading 2N in the Z-axis 
positive direction which act on TPM, the deformation effects 
of the liver model are shown in Fig. (3, b). In the MSM 
model, the deformation effects of the force point and edges 




are roughness, as shown in Fig. (3, a, b). But TPM model 
has better effects in deformation, especially the deformation 
of the force point and edges. 

2.2. Real-Time Simulation Condition 

We used the hardware platforms that included Omega. 7 
force-feedback devices which came from the Swiss Force 
Dimension Company, a computer which consisted of AMD 
Athlon 64X2 2.0GHz CPU 2.0GB RAM and Samsung LCD, 
as shown in Fig. (2, b). The graphic interfaces OpenGL can 
be used to display the effect of graphics rendering and visu- 
alization on this model simulation program. Programs of 
virtual scene setting and force-feedback effect came from the 
open source programs of the Swiss Force Dimension Com- 
pany (http://www.chai3d.org /do wnload.html). 

3. MODEL FORMULATION 

This model consisted of tension and three-parameter vis- 
coelastic structure [10] which had creep effects and viscoe- 
lastic properties. Finally, stress relaxation and creep formu- 
las had been obtained by mathematical formula derivation. 




oooooooo 



(a) (b) 
Fig. (1). (a) Layered sphere-filled structure, (b) The structure of layered sphere-filled model under the action of force. 




(a) 



(b) 




Fig. (2).(a)The structure of layered sphere-filled based on meshless method in TPM model has been established by using VC++ and the 
OpenGL. (b) Our test platform with a PC and a Omega. 7. 
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(a) 



(b) 



Fig. (3). (a) Loading 2N in the Z-axis positive direction, the liver model based on MSM presents the deformation rendering, (b) Loading 2N 
in the Z-axis positive direction, the liver model based on meshless method TPM model presents the deformation rendering. 
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Fig. (4). Three-parameter viscoelastic structure 



3.1. Model Framework 

A three-parameter viscoelastic structure consists of a 
Kelvin model and a spring as the basic model was used in 
this paper, as shown in Fig. (4). In Fig. (4), we assume that 

<J h (t) is the stress of three-parameter structure of two nodes 

at t time, E 1 and E 2 are the elastic modulus of the spring, 

7j l is the viscous modulus, € x (f) and £ 2 {t) are the length 
of the spring at t time. 

It will be damaged when the deformation length of the 
liver is more than 12mm [1 1] by using a rigid object to touch 
the liver surface with different velocity. So, the deformation 
length of TPM model should remain within 10mm as an ef- 
fective data by using a rigid sphere with a radius of 2mm to 
touch the model. Here, we set the radius of the sphere in 
TPM model as 10mm when exert external stress <j(t) is 

exerted to the outer sphere in this model, the stress status as 
shown in Fig. (5). It can be seen from Fig. (5) that the dis- 
placement of the outer sphere B changed under stress. To 
ensure the validity of the deformation, the displacement of 
sphere B should remain within 10mm. In order to simplify 
the calculation, we ignore the displacement of the spheres A 

and C. Here, we assume that A o is the center of the sphere 
A. Similarly, B o ,C o , D o are the centers of spheres B,C,D , 
respectively. And B o is the center of the sphere B after de- 



formation between two spheres which are connected by a 
three-parameter viscoelastic structure. 

We assume that s(t) is the distance between B o 

and B o . The initial length of line AB is x, and the length of 

line AB is y after deformation. The angle between x and y 
is 6 . Then, we assume As(t) = y — x . The length of defor- 
mation between A and B can be defined as 

0-cos# A 



As(t) = 8{t) 



sin# 



(1) 



The deformation length of this model should remain 
within 12mm to avoid the model damage. In order to ensure 

this criterion, 6 should satisfy that 6 e (0 — /r/6) . Here, we 
assume that cr f (t) is the tension. Then, the tension can be 
defined as 



a f (t) = 2Es(t) 



1 - cos 9 
sin# 



sin 6 



where E is the elastic modular. We 
sume/lj —2E (l-cos#), then Eq.(2) can be written 

a f {t) = X x -s{t). 

The formula of this model includes: 



(2) 

as- 

as 
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s(t) = s l (t) + s 2 (t) 

(T h (t) = E l £ l (t) + 7J l £ l (t) 
<J h (t) = E 2 £ 2 (t) 

a f (t) = A l -s(t) 

cr(t) = a h (t) + a f (t) 

By using inverse Laplace transform, we can obtain 
E x E 2 s(t) + E 2 ri x £(t) + \{E X + E 2 )s(t) + \rj x €{t) = 
(E x +E 2 )a(t) + ?] x &(t) 

This assumes 
7i 



(3) 
(4) 
(5) 
(6) 
(7) 

(8) 



P = 



p E l E 2 



E^ x l El 2 ^ 



E x +E 2 



can be written as 

(P 2 +A l )e(t) + (P 3 +\P l )m = <7(t) + P 1 <?(t) (9) 

Eq.(9) is the constitutive equation of this model. 

3.2. Relaxation 

After noting that s(t) = £ Q H(t) . Here, H(t) is the unit 
step function. It can be defined as 

f l, t>0 

H{t)- 

^0, t<0 

By using Laplace transform , Eq.(9) can be written as 



P 3 -P X P 2 , P x {P 2 +\) 



- + 



. + s 



(10) 



By using inverse Laplace transform, Eq.(10) can be writ- 
ten as 



cr(0 = (£ 2 +AK 



E 2 £ 



t{E x +E 2 )/ 



E x +E 2 



en) 



Eq.(l 1) is the relaxation equation of this model. 
3.3. Creep 

To exert a force O" 0 on soft tissue, there is a step- 
response function (7) . Considering loader^) = <J Q H(t) , 

substituting cf(t) = <y {) / S and &{t) = Sa(t) = <r () into 
Eq.(9), by using Laplace transform, Eq.(9) can be written as 

^ + P x cj, = {P 2 +\)s{s) + {P,+\P x )s8{s) (12) 
s 

Eq.(12) can be written as 



l + P x S 



S \P 2 +\) + {P i +\P l )S P 3 +\Py 



S(S+ ) 



+ - 



(p 2 +\) 



s + 



(13) 



This assumes 



P 2 +\ 1 



P 3 +M h 



— , by using inverse Laplace 



transform, Eq.(13) can be written as 



£(t) = - 



r l (l-e A ') + P l e' 



(14) 



Eq.(14) is the creep equation of this model. 
Eq.(14) shows that the sphere will be stable under the 
constant force O" 0 after t time, then the displacement of this 

sphere is s(t) . The relationship between O" 0 and s(t) is 

the stress-strain in this model. The displacement £ is 
changed with time t that can be defined as creep. Eq.(ll) 
shows that the force will change with the time when we give 
it a constant displacement £ Q . The force c(t) is changed 
with time t that can be defined as relaxation. 

4. RESULTS 

4.1. Model Parameters 

The effective elastic modulus of soft organ was estimated 
from the force versus displacement data of static indentations 
based on the small deformation assumption. Carter [12] and 
Kauer [13] used an instrument manually to measure the dis- 
placement and force response of the organ surface. In ex- 
periments with ex vivo pig liver, by using a hand-held probe 
Carter [12] calculated that the elastic modulus of pig liver 
which was equal to 490 kPa . Ottensmeyer [14] designed a 
robotic indenter to measure mechanical properties of pig 
liver during a minimally invasive surgery. The robotic in- 
denter can make small indentations in the range of 
±500 jUm . In experiments with in vivo pig liver, he calcu- 
lated that the elastic modulus of pig liver was equal to 2.2- 
15kPa. 

Effective shear and elastic moduli of pig liver were es- 
timated from the static indentation data by using the linear 
elastic contact theory and the small deformation assumption. 
The effective shear modulus, G, for small indentation of an 
elastic half space by a rigid hemispherical indenter can be 
obtained using the solution given by Lee [16]. So, by using 
the method of Lee the effective elastic modulus value was 
calculated between 11.6-18.7kPa [16]. On the other hand, 
in order to achieve better results, compared with the hand- 
held instrument, a robotic indenter was developed for mini- 
mally invasive measurement of the tissue properties. So we 
selected the parameter that was published by Ottensmeyer. 
Then we used the elastic modulus value of the liver which 
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was equal to 12 kPa in this paper by combining these pa- 
rameters. 

It can be seen from Fig. (5), under static stress <T 0 load- 
ing, we have s(t) = € x (f) + £ 2 (t) and 
<j Q = E x £ x (t) + E 2 £ 2 (t) . Then, we can write that 
cr 0 = Es(t) .Therefore, we get that E = E } = E 2 =12kPa. 

For A [ =2E(l- cos 6) , 0 e (0 - ;r/4) , the value of \ is 

calculated between 0.096kPa and 5.9kPa. 

Klatt [17] had concluded that the best agreement between 
fit and experiment was achieved with the Zener model. The 
reported value for viscous modulus base on Zener viscoelas- 
tic model was equal to 5.5 ± 1.6 Pa.s. Mohammad [18] cal- 
culated that the viscous modulus value was between 0.8-1.57 
Pa.s. Asbach [19] calculated that the viscous modulus value 
was equal to 5Pa.s. This deviation between our modeling 
results and other experiment results not only because of the 
model type (Voigt or Zener), but also because of different 
experimental conditions such as the slice thickness, kind of 
liver and in vivo or ex vivo execution. In this paper, the 
structure of TPM model was similar to that reported in As- 
bach [19], so we chose to use 5Pa.s as our viscous modulus 
value in TPM. 

4.2 Model simulation 

A force-feedback instrument was used to touch the 
No.71(x=0.002, y=0.136, z=0.057) sphere in the simulation 
experiments. The motion of other sphere was driven by the 
variation of the No. 71 sphere. And we selected above- 
mentioned parameters as the appropriate parameters for con- 
stitutive equations. 

Evren [15] used a probe with a diameter of 4mm to insert 
into the abdominal cavity of three different pigs. In the static 
loading test, the liver of each pig will indent slowly from 0 

jcr(t) 




to 8mm to characterize its force and displacement response. 
In Fig. (6), the blue curve, green curve and cyan curve repre- 
sent the force response value of the liver of pig #1, pig #2 
and pig #3, respectively. It can be seen from Fig. (6) that the 
force response value of three pigs for the indentation depth is 
different, but these data include a valid data range. In order 
to match the model experiment, a sphere with a diameter of 
4mm was used as a force-feedback instrument in TPM 
model, which is shown in Fig. (3, a, b). As shown in the Fig. 
(6), red curve represents the force response value of No. 71 
sphere in TPM for the indentation depth changing from 0 to 
8mm. The force response value of TPM model for the inden- 
tation depth maintains between the pig #1 and the pig #3. 

Noting that when t = 0 + in Eq. (1 1), we have 

a(0 + ) = (E 2 +A 1 )£ 0 (15) 

Here, we set that E 2 = \ 2kPa . 
Similarly, noting that when t — > go , we have 
E 2 £ 

a(co) = (E 2 +\)e 0 --^ r (16) 
tL x + E 2 

If the indentation depth of this model is 2 mm, as shown 
in the Fig. (5), the distance of two spheres is x, which will 
influence the model precision. When we assume that x is 

equal to 20mm, ^=0.119kPa can be obtained by using Eq. 

(2). Substituting \ into the Eq (11), the relaxation data in 

TPM model can be obtained. Similarly, when x is 25mm, 

30mm and 35mm respectively, \ will change as the value 

of x. It can be seen from Table 1, the result shows that both 
the model simulation data and the experimental data are very 
close. In the same way, when the indentation depth of this 
model is 4mm, 6mm and 8mm respectively, the relaxation 
data in TPM model is compared with the real liver experi- 
mental data, as shown in the Table 2, Table 3 and Table 4. 




Force response value of Pig liver #1 
Force response value of Pig liver #2 
Force response value of Pig liver #3 
Force response value of TPM 



0 1 2 3 4 5 6 7 

Displacement(mm) 



Fig. (5). The stress status of TPM model when exert the external 
stress to outer sphere. 



Fig. (6). The comparison between the force response of experimen- 
tal data of three pig livers and the force response of TPM model. 
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By analyzing these tables, we can find that the relaxation 
data in TPM are close to experimental results when 
x=30mm. Because the radius of the sphere in TPM is 
10mm, so we set that the distance of two spheres is 10mm. 

Here, we assume x=30mm and £ Q =4 mm, then the angle 

between two outer spheres should satisfy tan 0 = 4/30 , and 
we can get 6 = arctan(4 / 30) . Substituting 0 into the 
Eq.(2), /^=0.211kPa can be obtained. Finally, us- 
ing <j = F/Aq , we can obtain F=0.613N, which means the 
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instantaneous forces of liver. The relaxation curve is re- 
corded at the indentation depth of 4 mm at strain rates of 
4mm/s, 2mm/s, and 1 mm/s respectively in TPM, as shown 
in Fig. (10). 

Table 5 shows the comparison between the force re- 
sponse value of TPM and a valid data range of the force re- 
sponse value which was obtained from pig liver experiments. 
It can be seen from Tab. 5 that the force response value of 
TPM is in the range of the pig liver experimental data. So, it 
can be verified that the force response value in TPM was 
effective. 



Table 1. When the Indentation Depth is 2 mm, the Different Values of \ Were Obtained with Different x. The Table Shows the 
Comparison between The Relaxation Data of TPM Model and the Real Liver Experimental 



x (mm) 


\ (kPa) 


The Relaxation Data of TPM Within 30s (N) 


The Real Liver Experimental Relaxation Data Within 30s (AO 


20 


0.119 


0.304^0.153 


0.32^0.21 


25 


0.076 


0.303^0.152 


0.32^0.21 


30 


0.053 


0.302^0.152 


0.32^0.21 


35 


0.039 


0.301^0.151 


0.32^0.21 



Table 2. When the Indentation Depth is 4 mm, the Different Values of were Obtained with Different x. The Table Shows the 
Comparison between the Relaxation Data of TPM Model and the Real Liver Experimental 



x(mm) 


\ (kPa) 


The Relaxation Data of TPM Within 30s (N) 


The Real Liver Experimental Relaxation Data Within 30s (N) 


20 


0.466 


0.626^0.324 


0.62^0.32 


25 


0.301 


0.618^0.316 


0.62^0.32 


30 


0.211 


0.613^0.312 


0.62^0.32 


35 


0.155 


0.61^0.3 


0.62^0.32 



Table 3. When the Indentation Depth is 6 mm, the Different Values of \ were Obtained with Different x. The Table Shows the 
Comparison between the Relaxation Data of TPM Model and the Real Liver Experimental 



x(mm) 


\ (kPa) 


The Relaxation Data of TPM Within 30s (N) 


The Real Liver Experimental Relaxation Data Within 30s (N) 


20 


1.01 


0.98^0.528 


0.88^0.48 


25 


0.662 


0.95^0.5 


0.88^0.48 


30 


0.466 


0.93^0.48 


0.88^0.48 


35 


0.345 


0.928^0.47 


0.88^0.48 



Table 4. When the Indentation Depth is 8 mm, the Different Values of were Obtained with Different x. The Table Shows the 
Comparison between the Relaxation Data of TPM Model and the Real Liver Experimental 



x(mm) 


\ (kPa) 


The Relaxation Data of TPM Within 30s (N) 


The Real Liver Experimental Relaxation Data Within 30s(N) 


20 


1.71 


1.37^0.774 


1.72^0.68 


25 


1.14 


1.32^0.71 


1.72^0.68 


30 


0.81 


1.28^0.68 


1.72^0.68 


35 


0.6 


1.26^0.66 


1.72^0.68 
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Table 5. The Comparison Between the Force Response of TPM Model and a Valid Data Range of the Force Response was Ob- 
tained Through Pig Liver Experimental 



8 (mm) 


\ (kPa) 


The Value of Force in TPM Model (N) 


A Valid Data Range of the Force (AO 


2 


0.053 


0.152 


0.07-0.18 


4 


0.211 


0.312 


0.13-0.36 


6 


0.466 


0.48 


0.2-0.6 


8 


0.81 


0.68 


0.23-0.89 



In order to further verify the effectiveness of TPM, we 
compared the relaxation data in TPM with the relaxation data 
that was reported in Evren [15]. In stress relaxation experi- 
ments, the liver of each pig was indented to pre-defined 
depths of 2mm, 4mm, 5mm, 6mm, and 8mm respectively 
and recorded the force relaxation response value in the next 
30 seconds. In the Fig. (7), the blue curve, green curve and 
cyan curve represent the stress relaxation response value of 
pig #1, pig #2 and pig #3, respectively. The red curve repre- 
sents the stress relaxation response value of the #7 1 sphere in 
TPM at the indentation depth of 4mm with a strain rate of 
4mm/s in the next 30 seconds. It can be seen from Fig. (7) 
that the lines of stress relaxation value of TPM are close to 
the experimental data of pig liver. 

Here, we select that the force acting on the model is 0.1 96N. 
We assume /^=0.179kPa. Noting that when£ = 0 + , we 
have 



*(0 + ) = 



\+E 2 



(17) 



Then we can obtain the indentation depth of the model 
which is 1.31mm by substituting these parameters into 
Eq.(17). Similarly, noting that when£ — » oc , we have 



£(oo) = 



- + - 



(E x +E 2 )cj q 



(18) 



Then we can obtain the indentation depth of this model 
which is 3.8mm by substituting these parameters into 
Eq.(18). Similarly, when we select that the force acting on 

the model is 0.96N and we assume/^ =5. 8kPa by using 

Eq.(6). Noting that whenf = 0 + , we can obtain the indenta- 
tion depth of TPM which is 4.4mm by substituting these 
parameters into Eq.(17). Noting that when£ — » oo , we can 
obtain the indentation depth of TPM which is 10.86mm by 
substituting these parameters into Eq.(18). 

In order to further verify the effectiveness of TPM, we 
compared the creep data in TPM with the creep data that was 
reported in Amy [20] .In the Fig. (8), the cyan curve and the 
blue curve represent the lower boundary and upper boundary 
of creep data of the Ex vivo perfused pig liver under 0.1 96N 
respectively. The indentation depth of the Ex vivo perfused 
pig liver maintains between 3.8mm and 6.65mm. The red 
curve represents the creep data of the #71 sphere in TPM 
under 0.1 96N. It can be seen from the Fig. (8) that the creep 



curve in TPM is in the range of lower boundary to upper 
boundary. In the Fig. (9), the cyan curve and the blue curve 
represent the lower boundary and upper boundary of creep 
data of the Ex vivo perfused pig liver under 0.98N respec- 
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Fig. (7). The comparison between the force-time response of ex- 
perimental data of three pig livers and the force-time response of 
TPM model for the indentation depth of 4 mm with a strain rate of 
Amm/s in 30 seconds. 
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Fig. (8). The comparison between the creep response of Ex vivo 
perfused pig liver loading 0.1 96N in the Z-axis positive direction 
and the creep response of TPM model. 



A New Hybrid Viscoelastic Soft Tissue Model based on Meshless Method 



The Open Biomedical Engineering Journal, 2013, Volume 7 123 



20 



18 



0.98N Loading 

- Lower boundary of the creep response value of pig liver 
Upper boundary of the creep response value of pig liver 
■ Creep response value of TPM 



_ 14 




■ ooooooooooooooooooo 

^ 





w 5 


10 


15 


20 


25 


30 


35 



Time(s) 

Fig. (9). The comparison between the creep response of Ex vivo 
perfused pig liver loading 0.98Af in the Z-axis positive direction and 
the creep response of TPM model. 



4mm/s 
2mm/s 
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Fig. (10). The relaxation curve recorded at the indentation depth of 
4 mm at strain rates of 4, 2, and 1 mm/s in TPM model. 

tively. The indentation depth of the Ex vivo perfused pig 
liver maintains between 9.6mm and 11.9mm. The red curve 
represents the creep data of the #71 sphere in TPM under 
0.98N. It can be seen from Fig. (9) that the creep curve in 
TPM is in the range of lower boundary to upper boundary. 

5. DISCUSSION AND CONCLUSIONS 

Comparing with the experimental data of pig liver which 
was reported in Evren [15] and Amy [20], we verified that 
the TPM was effective in stress-strain, stress relaxation and 
creep data. Moreover, the TPM with two-layer spheres struc- 
ture had been established. On the basis of this model, in or- 
der to further improve the model precision, we can establish 
TPM structure with three or several layer spheres. 

Although this model had stress-strain, stress relaxation 
and creep deformation properties, which had a significant 
advantage over the MSM, but the precision of this model 



needs to be improved further. And the data gained from the 
simulation experiment had some deviation. In the future, we 
will refine the model to increase the precision of the model, 
and the model will be evaluated on more challenging and 
complicated tasks. 
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